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Abstract. We study the dependence of the quantum graph HamHtonian, its resolvent, and its spec- 
trum on the vertex conditions and edge lengths. In particular, several results on the analyticity and 
interlacing of the spectra of graphs with different vertex conditions are obtained and their applications 
are discussed. 
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■ 1. Introduction 

Q Graph models have long been used as a simpler setting to study complicated phenomena. Quantum 
^ graphs in particular have recently gained popularity as models for thin wires, eigenvalue statistics 
I — I of chaotic systems and properties of the nodal domains of eigenfunctions. We refer the interested 
^ ^ reader to the recent reviews [ini[I3,ll9] and collections of papers [3l|7j. 

A quantum graph is a metric graph equipped with a self-adjoint differential "Hamiltonian" operator 
(usually of Schrodinger type) defined on the edges and matching conditions specified at the vertices. 
Every edge of the graph has a length assigned to it. In this manuscript we establish several results 
, G^ . concerning the general properties of the spectrum of the Hamiltonian as a function of the parameters 
'. involved: the edge lengths and matching (vertex) conditions. In Section [2] we introduce the standard 
^ ■ notions related to quantum graphs, vertex conditions, and quadratic forms of the corresponding 
Hamiltonians. The most widely used example of the vertex conditions, the 5-type condition, is 
)0 . described in some detail. The results presented in Section [3] concern the analyticity of the spectrum 
as a function of the parameters involved: vertex conditions and edge lengths. Analytic dependence 
on the potential (as an infinite-dimensional parameter) can also be established by the methods used 
O ■ in the Section, but we omit this to keep the presentation simple. In Section |4] we compute the 
O . derivative of an eigenvalue with respect to a variation in the length of an edge; this is an analogue 
. . of a well-known Hadamard variational formula. Derivative with respect to a parameter of a 5-type 

[ condition is also computed. 
^ Section [5] again focuses on the (5-type condition at a vertex. Varying such a condition at a specified 

■ vertex or changing connectivity of the vertex we obtain families of spectra and study their interlacing 
properties. Eigenvalue interlacing (or bracketing) is a powerful tool in spectral theory with such well- 
known applications as the derivation of the asymptotic Weyl law, see [5]. In the graph setting, it 
allows one to estimate eigenvalue of a given graph via the eigenvalues of its subgraphs, which may 
be easier to calculate. Interlacing results on graphs have already been used in several situations 
[2l|2Tl|26]. We significantly generalize these results and put them in the form particularly suited for 
the applications. We discuss several applications. In particular, we give a simple derivation of the 
number of nodal domains of the n-th eigenfunction on a quantum tree and study irreducibility of the 
spectrum of a family of graphs obtained by varying 5-tjpe condition at one of the vertices. 

2. Quantum graph Hamiltonian 

Let F = (V, S) be a graph with finite sets of vertices V = {vj} and edges S = {ej}. We will assume 
that r is a metric graph, i.e. every edge e is a 1-dimensional segment with a positive finite length 
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Lg and a coordinate assigned. Each edge corresponds to two directed edges, or bonds, of opposite 
directions. If bonds hi and 62 correspond to the same edge, they are called reversals of each other. 
In this case we use the notations hi = 62 and 62 = A bond h inherits its length from the edge e 
it corresponds to, in particular, = = L^. A coordinate Xf, is assigned on each bond, and the 
coordinates on mutually reversed bonds are connected via xj; = Lf, — Xb- 

Definition 2.1. 

• The space L2{T) on F consists of functions that are measurable and square integrable on each 
edge e with the norm 

ii/iiW) :=i:ii/iii(e)- 

In other words, L2{T) is the orthogonal direct sum of spaces L2{e). 

• We denote by H''(T) the space 

^^(r):=0i7^(e), 

which consists of the functions / on F that on each edge e belong to the Sobolev space H^{e) 
and is equipped with the norm 

ll/ll|rfc(r) X/ ll.^ll?/'=(e)- 

• The Sobolev space H^iT) consists of all continuous functions from H^{r). 

Note that in the definition of if^(F) the smoothness is enforced along edges only, without any 
junction conditions at the vertices at all. The continuity condition imposed on functions from the 
Sobolev space H^{T) means that any function / from this space assumes the same value at a vertex 
V on all edges adjacent to v, and thus f{v) is uniquely defined. This is a natural condition for one- 
dimensional if^-functions, which are known to be continuous in the standard 1-dimensional setting. 

A metric graph becomes quantum after being equipped with an additional structure: assignment of 
a self-adjoint differential operator. This operator will be also called the Hamiltonian. The frequently 
arising in the quantum graph studies operator is the negative second derivative acting on each edge 
{x is the coordinate x along an edge) 



(1) f{x) ^ 

or the more general Schrodinger operator 



dx"^ 



(2) + 

where V{x) is an electric potential. (Our results will hold, for instance, for V e L2{T). Generaliza- 
tions to more general operators, e.g. including magnetic terms, are also straightforward.) 

Notice that for both these operators the direction of the edge is irrelevant. This is not true anymore 
if one wants to include derivative term of an odd order, e.g. magnetic potential, but we shall not 
address such operators in the present note (see, e.g., [8l|2l] concerning these issues). 

The natural smoothness requirement coming from the ODE theory is that / belongs to the Sobolev 
space -ff^(e) on each edge e. Appropriate boundary value conditions at the vertices {vertex conditions) 
still need to be added, which are considered in the next subsection. 
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2.1. Vertex conditions. We will briefly describe now the known descriptions of the vertex condi- 
tions that one can add to the differential expression ([2]) in order to create a self-adjoint operator (see, 
e.g., [H1II31[T61[I8] for details). 

Assume that the domain of the operator is a subspace of the Sobolev space H'^(T) (see the references 
above for the justification of this assumption). Then the standard Sobolev trace theorem (e.g., |6]) 
implies that both the function /, / G -f^^(e), and its first derivative have correctly defined values at 
the endpoints of the edge e. Thus, for a function / G H'^^T) and a vertex v we can define the column 
vectors F{v) and F'{v) 



(3) 



F{v) :- 



F'iv) :-- 



of the values at the vertex v that functions / and /' attain along the edges incident to v. Here is 
the degree of the vertex v. The derivatives of / at vertices are always taken away from the vertices 
and into the edges. 

Descriptions of all possible vertex conditions that would make operator l-L self-adjoint can be done 
in several somewhat different ways. Below we list the most usual descriptions, which were introduced 

in [I3lESl[IB]. 

Theorem 2.2. Let T he a metric graph with finitely many edges. Consider the operator l-i acting 

as + on each edge e, with the domain consisting of functions that helong to H'^{e) and 

satisfying some local vertex conditions involving vertex values of functions and their derivatives. The 
operator is self-adjoint if and only if the vertex conditions can he written in one (and thus any) of 
the following three forms: 

A: For every vertex v of degree there exist x d^ matrices and By such that 



(4) 



The d^ X 2d^ matrix {A^ B^) has the maximal rank. 



(5) The matrix A^B* is self-adjoint. 

and functions f from the domain of l-L satisfy the vertex conditions 

(6) AyF{v) + ByF'{v) = 0. 

B: For every vertex v of degree d^, there exists a unitary dy x d^ matrix Uy such that functions 
f from the domain of Ti satisfy the vertex conditions 

(7) i{Uy - I)F{V) + {Uy + I)F\v) = 0, 

where I is the dy x dy identity matrix. 
C: For every vertex v of degree dy, there are three orthogonal (and mutually orthogonal) projec- 
tors Pd,v, Pn,v and Pji^y := I — Pd,v — Pn,v (one or two projectors can he zero) acting in C^" 
and an invertihle self-adjoint operator Ay acting in the subspace Pr^^C^", such that functions 
f from the domain of l-L satisfy the vertex conditions 

{Pj^ yF{v) = - the "Dirichlet part", 
Pn,vF'{v) = - the "Neumann part", 
Pji^yF'{v) = KyPR.,F{v) - the "Rohm part". 
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2.2. Quadratic form. To describe tlie quadratic form of the operator "H, which is a self-adjoint 
reahzation of the Schrodinger operator ([2]) acting along each edge, the self-adjoint vertex conditions 
written in the form (C) of Theorem 12.21 are the most convenient. The following theorem is cited 
from [TSl. 



Theorem 2.3. The quadratic form h ofH is given as 



(9) 



where (, ) denotes the standard Hermitian inner product in C^™^^-" . The domain of this form consists 
of all functions f that belong to H^{e) on each edge e and satisfy at each vertex v the condition 
Pd,vF = 0. 

Correspondingly, the sesqui-linear form ofH is 



(10) 



h[f,9] 



ee£ 



I ^%dx + Y,I V{x)f{x)g{x)dx + J2 {^vPr,.F, Pr,,G) . 



2.3. Examples of vertex conditions. In this paper we will often be dealing with the 6-type 
conditions which are defined at a vertex v as follows: 

f{x) is continuous at v, 

where for each vertex v, is a fixed number. One recognizes this condition as being an analog of 
the conditions one obtains for the Schrodinger operator on the line with a 6 potential. The special 
case = is known as the Neumann (or Kirchhoff) condition. 
The 5-type condition can be written in the form with 



/ 1 




A,, 





-1 











a. 



and 



/O 



B., 



1 


0\ 



-1 
0/ 
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the self-adjointness condition ([5]) is satisfied if and only if is real. 

In order to write the vertex conditions in the form ([8]), one introduces the orthogonal projection 
Pd,v onto the kernel of B^, the projector Pr^^, = I — Po,v and the self-adjoint operator A^, = B^ ^^A^ 
on the range of Pr^v A straightforward calculation shows that Pr,„ is the one- dimensional orthogonal 
projector onto the space of vectors with equal coordinates and thus the range of Pd,v is spanned by 
the vectors rk, k = 1, ...,dy — l, where has 1 as the fc-th component, —1 as the next one, and zeros 
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a, 



otherwise. Then A^, becomes the multiphcation by the number — . In particular, the quadratic form 

dv 

of the operator H (assuming 5-type conditions on all vertices of the graph) is 



h[f, f] 



(12) 



E 



df 



dx 



df 



dx 



dx + J2j V{x)\f{x)fdx + J2{L,F,F) 
Vix)\fix)\'dx + J2^v\fiv)\' 



dx + 



defined on / G H^(T), which are automatically continuous, and so F{v) = (/(f), f{v)Y. 

Vertex Dirichlet condition requires that the function vanishes at the vertex: f[v) = 0. At the 
first glance, it might look like it is significantly different from the 5-type conditions, but a closer 
inspection shows that this is not the case. Indeed, since the function must vanish when approaching 
the vertex from any edge, the vertex Dirichlet condition can be recast in the following form: 



(13) 



f{x) is continuous at v, 
f{v) = 0, 



Now one finds resemblance with ( II ip . and indeed, if one divides the equality in (ITT]) by and then 
takes the limit when — >■ oo, one arrives to (I13p . 

Hence, vertex Dirichlet condition seems to be the limit case of ( II ip when — )■ oo. We thus 
introduce the extended 6-type conditions by allowing = oo. In order to avoid considering infinite 
values of the two types of conditions can be also written in the form 



(14) 



cos 



Here 7t, = corresponds to the Neumann condition and 7^, = 7r/2 corresponds to the Dirichlet one, 
with more general (5-type conditions in between. Usefulness of considering Dirichlet condition as a 
part of the family of 5-type conditions becomes clear in spectral theory, as will be illustrated, for 
instance, in Theorem 15.11 

In fact, it will be convenient later on in this paper to rewrite the extended 5-type conditions ([E 
in the following form: 



(15) 



1. 



where z belongs to the unit circle in the complex plane, i.e. \z 

Interpreting the Dirichlet condition in terms of the corresponding projectors, as in part C of 

I and, correspondingly, Pr, 



Theorem 12. 2[ one notices that here Pd,v = I and, correspondingly, P^^, = 0. Hence there is no 
additive contribution to the quadratic form h[f, f] coming from the vertex v. Instead, the condition 
f{v) = is introduced directly into the domain D{h). 

To summarize, the quadratic form for a graph with the extended 5-type conditions (i.e., allowing 
a,, = 00) at all vertices, can be written as 



(16) 



where 



df 



dx 



dx 



E 



V{x)\f{x)\^dx+ J2 

{veV I a„<oo} 



/ G H^(T), f is continuous on F, and f{v) = whenever = 00. 
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Vertex Dirichlet condition is an example of a decoupling condition, since it essentially removes any 
connection between the edges attached to the vertex. Another example that will be useful to us is 
Pr,v = I, Pn,v = Pd,v = and = diag(ai, . . . , a^^). In this case the function is no longer required 
to be continuous at the vertex and the condition reduces to 

on every edge e incident to the vertex v. 

3. Dependence on vertex conditions and edge lengths 

3.1. Dependence of the Hamiltonian. In this section, we discuss the (analytic) dependence of the 
quantum graph Hamiltonian "H on the vertex conditions and the edge lengths. This issue happens 
to be important in many circumstances, e.g. when considering dependence of the spectrum and 
the eigenf unctions. To keep the notation simpler, we only address the case V = 0. However, the 
results can be extended to non-zero electric potentials without any change in the proofs. In fact, the 
dependence on the potential is also analytic, so the potential can be added to the vertex conditions 
and edge length as an extra (infinite-dimensional) parameter. 

As before, the graph F is assumed to be finite (i.e., it has finitely many vertices and finite lengths 
edges) . 

It will be convenient here to consider the vertex conditions in the form of (Q: 

A,F{v) + B.F'iv) = 0. 

Given a set of matrices = {A^ By)\y^\;, we denote by A the collection of Ay for all v eV: 

A:= vev- 

Then A can be considered as a block-diagonal matrix of the size {J2v '^f) ^ '^iJ2v '^f )' with individual 
blocks of sizes x 2dv. The space of such complex matrices A can be identified with C^^^'^^, or 
just C^"^, where we will use the shorthand notation 

V 

We now define the sub-set U of C^'^ that consists of matrices A satisfying the maximal rank condition 
dlD for each f e V: 

U := {A I {A^ B^) has maximal rank for any v G V}. 
We denote by Us the subset of U consisting of matrices satisfying the self-adjointness condition (E}: 

Us:= {AeU \ A.^B* is self-adjoint for any v e V}. 
The following statement describes some simple properties of these sets: 

Lemma 3.1. 

(1) The complement C^^ \ UofUin C^'' is algebraic. 

(2) U is an open and everywhere dense domain of holomorphy in C^*^. 

Proof. Algebraicity of the set C^'' \ U is clear, since its elements are described by the algebraic 
relations forcing the highest order minors to vanish. This proves the first statement of the Lemma. 
The second claim is an immediate corollary of the first one, if one can show that U is not empty. 
This is done by noticing that A := (I 0) E U. □ 

Let us now return to the quantum graph Hamiltonian. Since we are going to look into its depen- 
dence on vertex conditions, we introduce the corresponding notation: 
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Definition 3.2. The operator —(P/dx^ defined on functions from H'^iT) satisfying ([6]) at eacli vertex 
V is denoted by "H^^. 

As we liave already seen, tlie domain V{T-L^) of tlie operator Ti^ is a closed subspace of H'^{T) 
tliat can be described as follows: 

:= {/ G H\T) I A^F{v) + B^F'iy) = for all v G V}. 

In other words, "D^ is the kernel of the continuous linear operator 

Ta : ^'(r) ^ C'^, 

where 

(17) TA-.f^ {A,F{v) + B^F'iv)} G C^" = C^^. 

This simple observation allows us to establish nice behavior of the domain of "H^ with respect to the 
vertex conditions matrix A. In order to do this, let us start with a simple lemma: 

Lemma 3.3. 

(1) The operator function A ^ Ta is analytic in C^'^ with values in the space L{H^(T),C^^) of 
bounded linear operators from H'^{T) to C^"^. 

(2) For any A E U, the operator Ta is surjective. 

Proof. Indeed, according to ( IT71) . the function is in fact linear, and thus analytic with respect to A. 
Also, the set of vectors {F{v), -F'('y)}|i,gv achievable from elements / of ©e-f^^(e) is clearly arbitrary. 
Then, if the maximal rank condition is satisfied, this implies the surjectivity of Ta- CH 

Let us consider the trivial vector bundle 

U X H^{T) U 

over U with fibers equal to if^(r). Consider the sub-set 

V :=[j {{A} xVa)(1Ux H\V). 

In other words, we look at the domain Va of the operator TIa as a "rotating" with A subspace 
of if^(r). The next results shows that this domain rotates "nicely" (analytically) with A and is 
topologically and analytically trivial as a vector bundle. 

Theorem 3.4. 

(1) V is an analytic sub-bundle of co- dimension 2E of the ambient trivial bundle. 

(2) The bundle V is trivializable. In other words, there exists a trivialization, i.e. an analytic 
operator- function T{A) on U with values in linear bounded operators from a Hilbert space H 
into H'^(T), such that kerT(^) = and ranT(^) = Va for any A E U. 

(3) After the trivialization, the values of the resulting analytic in U operator- function 

A ^ HaT{A) 
are Fredholm operators of index zero from H to L?'{T). 

Proof. The first statement of the Theorem is local. So, let us pick a matrix ^ U . According to 
Lemma [3l3l the operator Tao is surjective, and thus has a continuous right inverse R. Thus, T^q-R = I. 
Then TaR is invertible for A close to ^o- This implies that for such A, one has TaR{TaR)~^ = I- 
This implies that 

P{A):=I-R{TaR)-'Ta 
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is a projector onto the kernel of 7^, i.e. on V^x- Since this projector, by construction, is analytic 
with respect to ^ in a neighborhood of Ao, this proves a part of the first claim of the theorem: 1) is 
an analytic Hilbert sub-bundle (e.g., [2Z])- It only remains to notice that co-dimension of the kernel 
of 7^ is the dimension of its range, i.e. 2E. 

To prove the second claim, we notice that V is an infinite dimensional analytic Hilbert bundle. 
Due to the Kuiper's theorem on contractibility of the general linear group of any infinite-dimensional 
Hilbert space [20], all such bundles are topologically trivial. Since the base U is holomorphically 
convex, the Bungart's theorem |3] says that the same holds in the analytic category (see further 
discussion of the technique and relevant references in the survey [27]). 

Let us prove the third claim. Since Ha is the restriction to V^, of the fixed operator —(P/dx^ (with 
no vertex conditions attached), acting continuously from H^{T) to L'^iV), the operator-function in 
question can be written as {—(P/dx'^)T{A) and thus is analytic. The Fredholm property and the 
zero value of the index follow from [SI Theorem 14]. □ 

3.2. Dependence of the resolvent. We now consider the question about the dependence on vertex 
conditions of the resolvent of H^. In order to do so, we need to consider the operator family — il, 
where I denotes the identity operator in L'^iT)- Sobolev's compactness of embedding theorem shows 
that I is a compact operator from H'^{T) to L2{T). This and the previous theorem imply that "H^ — il 
is also an analytic family of Fredholm operators of zero index. 

As analytic Fredholm theorem (see Appendix shows, the set S of the matrices A for which the 
operator "H^ — il : — ^ -^^(r) is not continuously invertible, is principal analytic (i.e., can be given 
by an equation 0(^) = 0, where function is analytic). Since this set S does not include any points 
A of Us (because in this case the operator is self-adjoint), this singular set is nowhere dense. 

Theorem 3.5. The resolvent ("H^ — if)^^ , defined in U \ T,, is analytic and has values that are 
compact operators in L'^(T). 

Proof. Analyticity of the inverse to an analytic family of bounded invertible operators is well known 
(e.g., [I5l|27]). Thus, {Ha ~ "^I)"^ is an analytic family of operators from L'^{T) to iJ^(r), and thus 
also as a family of operators acting in L'^{T). However, considered as operators in L^, they factor 
through the compact embedding of into L'^{T), and thus are compact. □ 

3.3. Variations in the edge lengths. Sometimes one needs to consider the quantum graph's de- 
pendence on the variations in the edge lengths parameters {^e} (without changing graph's topology). 
We thus extend the previous considerations to include the dependence on the vector 

in fact, we allow "complex values of lengths" £, i.e. £ G (C \ {0})-^. 

Let r be a quantum graph with the Hamiltonian = —d'^/dx^ equipped with vertex conditions 
described by a matrix A = {{A^ B^)\y^\;}. We would like to vary the lengths of the edges (indepen- 
dently from each other), without changing topology or vertex conditions. Let us consider the vector 
^ = {^e} G (M"'")'^ of dilation factors along each edge. It is clear that such a dilation is equivalent to 
keeping the metric graph structure the same, while replacing "H = —d'^/dx'^ with "H^ = {— ^"^rf^/rfx^ 
and {Ay By) with {Ay BySy), where is the diagonal dy x dy matrix having C,e as its diagonal entries, 
where e denotes the edges incident to the vertex v. 

Definition 3.6. For any A from the previously described set U of complex matrices A = {{Ay By)} 
satisfying the maximal rank condition and for any vector ^ G (C\{0})^, we denote by 'H^,^ the Hamil- 
tonian —^~'^d'^/dx1 on r with the vertex conditions provided by the matrices ^4.^ := {{Ay By'Ey)}. 
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Notice that this rescahng allows for "complex lengths" of the edges. This is useful, when considering 
analytic properties of the Hamiltonian with respect to the parameters. 

The same arguments as in proving Theorem 13.41 lead to the following statement: 

Theorem 3.7. 

(1) The bundle V over U x {C \ {0})^ with the fiber over (A,^) defined by the vertex condition 
matrix A^, is an analytic sub-bundle of co-dimension 2E of the ambient trivial bundle 

(f/ X (C \ {0})^) X H\T) ^ (f/ X (C \ {0})^) . 

(2) The bundle T> is trivializable. In other words, there exists a trivialization, i.e. an analytic 
operator-function T{A,S,) on U x [C\ {O})'^ with values in linear bounded operators from 
a Hilbert space H into H'^{r), such that kerT(^, ^) = and ranT(^, = for any 

(AOet/x(C\{0})^. 

(3) After the trivialization, the values of the resulting analytic in U x (C\{0})'^ operator-function 

A^'Ha,^T{A,0 
are Fredholm operators of index zero from H to L^{T). 

3.4. Dependence of the spectrum on the vertex conditions. In this section, we will derive 
some basic facts about relations between the quantum graph eigenvalues and continuous graph pa- 
rameters (such as matrices A = {(^^ 6i,)}|«ev of vertex conditions and edges' lengths {Le}\ei^s- We 
denote by 'H_4 the operator —d'^/dx'^ on a compact graph F with the domain described by the vertex 
conditions ^ that correspond to the matrix A = {{A^ bv)}\vev C U. Here U C C^'' (with q = J2^v) 
consists of such complex A that the rank of the x 2(i^-matrix {A^ B^) is maximal for any vertex 
V. We will also use notation 'Ha,£_ for the re-scaled version of "H^ that acts as —^'"^d^/dx^ on the 
domain described by A^^ = (A^, B^E^). Here ^ = {^e} ^ (C \ {0})^ is the vector of (non-zero) scaling 
factors applied on each edge. 

We are interested in the dependence of the spectrum of TIa,^ on the parameters {A,C,) C U x {C\ 
{0})^. Thus, we consider the operator pencil 71^,^ — AI of unbounded operators in L'^iT). 

The following result, which addresses analytic behavior of the spectrum, follows from Theorem 13.71 
and the analytic Fredholm theorem (Appendix R|) and [27, Theorem 4.11]. 

Theorem 3.8. 

(1) The set S of all vectors 

{A,^,X) cUx (C\{0})^ X C 

such that Ha,^ — AI does not have a bounded inverse in L'^(T), is principal analytic. Namely, 
there exists a non- zero function ^{A,^, A) analytic in Ux (C\{0})'^xC, such thatS coincides 
with the set of its zeros. 

(2) For any integer k > 1, the set Sk of all vectors 

{A,^,\) cUx (C\{0})^ X C 
such that dim Ker {'Ha,i^ — AI) > k, is analytic. 
Remark 3.9. 

• The set S is the graph of the multiple-valued function 

(AO^^(^AC), 

and can be considered as a kind of "dispersion relation." Here ct('H) denotes the spectrum of 
the operator T-i. 
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• The sets Sk are clearly nested: iS^+i C Sk- Moreover, Si = S. 

• For A E Us (i.e., satisfying the condition guaranteeing self-adjointness of H^) and ^ G 
(M \ {0})^, the operator 'Ha,£, is selfadjoint and its spectrum is real and discrete. 

3.5. Eigenfunction dependence. We present here, for completeness, a simple and well known 
consequence of perturbation theory. 

Theorem 3.10. The kernel Ker{'Hji^^^ — AI) forms a holomorphic k- dimensional vector bundle over 
the set Sk \ Sk+i- In other words, if one has an local analytic eigenvalue branch A(^, ^) of constant 
multiplicity k, then, at least locally, one can choose an analytic basis of k eigenf unctions. 

Proof. The spectral projector onto the corresponding spectral subspace is clearly analytic with respect 
to the parameter. □ 

4. An Hadamard type formula 

Hadamard's variational formulas |9l[l2l[Tll|22] deal with the variation of the spectral data with 
respect to the domain perturbation. For simplicity, we will consider the case of variation of the 
eigenvalues with respect to a change in a "loose" edge's length. Namely, the end of the edge is 
assumed to be a vertex of degree 1 and the length of the edge will be denoted by s. We will also use 
s to denote the vertex of degree 1. The proof follows the well established pattern (see, e.g. [H]) and 
can be easily generalized. 

Proposition 4.1. Let A = A(s) be a simple eigenvalue of a graph F with a loose edge of length s and 
the Dirichlet condition imposed at the end-vertex s. Let f = fs{x) be the corresponding normalized 
eigenfunction. Then 

(18) S-I^'WI^ 

where f'{s) is the value of the derivative of the eigenfunction at the end-vertex. 

Proof. First we remark that by preceding theorems we can differentiate both the eigenvalue and the 
eigenfunction. We will omit the subscript s of f unless we want to highlight the dependence of the 
eigenfunction on s. 

The Dirichlet condition at the vertex s has the form 

fsis) = 0. 

Differentiating it with respect to s we get 

(19) ^ + f{s)=0. 

On the other hand, the eigenfunction is L2-normalized. Differentiating the normalization condition 
(/,/) = Owe get 

2<Ke(/.f)=0. 

where we used the fact that the loose edge's contribution to the derivative is 

(20) \fs{x)\'dx = \f{s)f + 2Dltl^ f^dx 

and /s(s) = 0. Finally, representing the eigenvalue as A = h[f, /], where h is the quadratic form (|9]), 
we get 

Of 



(21) ^ = \f^s)f + V{s)\f{s)f + 2mch 



f 1 o 

OS 
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To evaluate the sesqui-linear form h we note that the derivative df/ds satisfies the same vertex 
conditions as / everywhere apart from the point s. Integrating by parts we get 



h 



ds 



dl 

ds 



Now we use that "H/ = A/ and equations ( |T9|) and ( l20l) to get 



2^th 



f: 



df 
ds 



ds 



2\f'is)\' = -2\f'is) 



Substituting the last equation and the Dirichlet condition into fj2T|) yields the desired result. □ 

One can obtain similar results when varying vertex conditions rather than the length. For simplicity 
we only consider the S-tjpe vertex conditions. 

Proposition 4.2. Let A = X{a) be a simple eigenvalue of a graph T which satisfies S-type vertex 
condition at v with the parameter a ^ oo. Then 

(22) ^ = \m\'. 



da 



If we re-parameterize the conditions at v as 
(23) 



-fiv), 



now allowing Dirichlet (( = 0) and excluding Neumann (( = oo) conditions, the derivative is 

2 

(24) 



dX 
d( 



Proof. The proof follows the pattern of the proof of Proposition 14.11 with minor modifications. We 
deal with a-derivative first. The derivative of the normalization condition is now 



0. 



Taking the derivative of the quadratic form, see (fT2l) . yields 



^ = \fiv)\' + 2D\th 
da 



J } rN 

oa 



Integrating by parts inside the sesqui-linear form h and collecting together all the "boundary" terms 
at f , we get 



h 



f ^ 
oa 



«l^(^)/(^) 



^(^) y ^(^) + (nf ^) = (uf ^ 

da ^ dxe V ^ da ) \ 'da 



where we used the S-tjpe condition at the vertex v. Now we use Tif = \f and the derivative of the 
normalization condition to get 

2D\th /,|^ =0, 
da 



and, therefore, 
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To deal with the Dirichlet case, we calculate, for a 7^ and ( ^ 0, 



2 



using condition ( !23l) in the last step. Since X{() is an analytic function, the value of the derivative 
at C = now follows by continuity. □ 

5. Eigenvalue interlacing 

We will be assuming that the Hamiltonian T-L is —d'^/dx'^ + V{x) with the vertex conditions specified 
in the results. The eigenvalues A„ of T-L (also referred to as the eigenvalues of the graph and denoted 
A„(r)) are labeled in non- decreasing order, counting their multiplicity. 

The first theorem of this section describes the effect of modifying the vertex condition at a single 
vertex. We denote by a compact (not necessarily connected) quantum graph with a distinguished 
vertex v. Arbitrary self-adjoint conditions are fixed at all vertices other than v, while v is endowed 
with the 5-type condition with coefficient a: 

f is continuous at v and 

e€£v 

Theorem 5.1. Let F^' be the graph obtained from the graph by changing the coefficient of the 
condition at vertex v from a to a' . //— 00 < a < a' < 00 (where a' = 00 corresponds to the Dirichlet 
condition, see section \279\] . then 

(25) A„(F,) < A„(F,0 < A„+i(F,). 

// the eigenvalue A„(Fq,') is simple and it's eigenfunction f is such that either f{v) or ^ /'(f ) is 
non-zero then the inequalities can be made strict, 

(26) A„(F,) < A„(F,0 < A„+i(F,). 

Proof. The case of strict inequalities follows simply from the positivity of the derivative of A„(Fa') 
with respect to the parameter of the vertex condition. Proposition 14.21 For the possibly degenerate 
case we directly use the monotonicity of the quadratic form and rank-one nature of the perturbation. 

Denoting by Fgo the graph with the Dirichlet condition at the vertex v, we will actually prove the 
chain of inequalities 

(27) A„(F„) < A„(F,0 < A„(Foo) < A„+i(F,), 

which is obviously equivalent to inequality (1251) . Since we are now considering the Dirichlet case 
separately, we will assume that a' ^ 00. 

Consider the quadratic forms ha, ha' and h^o of the corresponding Hamiltonians. According to 
the discussion in section 12. 3[ we have 

ha[fj]=h^[fj]+a\fiv)\' and ha'[fj] = h^[fj] + a'\fiv)\' 

on the appropriate subspaces of H^{T){= 0egf -f^^(e)). In fact, D{ha) = D{ha') and D{hoo) = {/ G 
D{ha) : f{v) = 0}. 

All inequalities follow from the min-max description of the eigenvalues, namely 

(28) A„ = min max h\f, f]. 

XcD:dim(X)=n /gX: ||/|| = 1 

The first inequality in (|27|1 follows immediately from the observation that ha' > ha for all /. 
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The domain D{hoo) is smaller than D{ha') and the forms hoo and ha' agree on D{h^). Minimization 
over a smaller space results in a larger result, implying the second inequality in ( l27l) . 

The last inequality follows from the fact that D{h^) is a co-dimension one subspace of D{ha')- 
To provide more detail, let the minimum for Xn+i{Ta) be achieved on the subspace Y (which is the 
span of the first n + 1 eigenvectors) of dimension n + 1. Then there is a subspace Y^o of dimension 
n, such that Y^o C Y and Y^ C D{h^). Then 

An(roo) = niin max/ioo < max /iqo 

X:dim(X)=n fdX f&oo 

= max ha < max ha = \n+i(Ta)- 
This is precisely the last needed inequality. □ 



This theorem allows us to prove a simple but useful criterion for the simplicity of the spectrum of 
a tree. 

Corollary 5.2. Let T he a tree with a S-type condition at every internal vertex and an extended 
6-type condition at every vertex of degree 1. If the eigenvalue X of T has an eigenfunction that is 
non-zero on all internal vertices ofT, then A is simple. 

Equivalently, if an eigenvalue X of the tree T is multiple, there is an internal vertex v such that all 
functions from the eigenspace of X vanish on v. 

Proof. The two statements are almost contrapositives of each other, modulo the following observation: 
if for every internal vertex there is an eigenfunction that is non-zero on it, one can construct an 
eigenfunction which is non-zero on all internal vertices. Indeed, if m is the dimension of the eigenspace 
of A, then the subspace of the eigenfunctions vanishing on any given v is at most m — 1. A finite 
union of subspaces of dimension m — 1 cannot cover the entire eigenspace. 

We will work by induction on the number of internal vertices. If a tree has no internal vertices, it 
is an interval and there is nothing to prove since all eigenvalues are simple. 

Assume the contrary: there is an eigenfunction / which is not zero on all internal vertices of the 
tree, but A is not simple. Take an arbitrary internal vertex v and another eigenfunction g. Cutting 
the tree at the vertex v we obtain dy subtrees. On at least one of them the function g is not identically 
zero and not a multiple of /. Let T' be a such a subtree. 

Then there is an a < oo such that (A, /) is an eigenpair of the tree T^^ endowed with the condition 
f'{v) = af{v). Similarly, there is an a' < cxd such that (A, (7) is an eigenpair of the tree T'^i. By 
inductive hypothesis, A is simple on T^, therefore a' 7^ a. This, however, contradicts inequality (|26|) . 

□ 



The next theorem deals with the modification of the structure of the graph by gluing a pair of 
vertices together. 

Theorem 5.3. Let T he a compact (not necessarily connected) graph. Let vq and vi he vertices of 
the graph T endowed with the 6-type conditions, i.e. 



f is continuous at Vj and 

E i(^.) = «.M), J = 0,1. 



Arbitrary self-adjoint conditions are allowed at all other vertices ofT. 
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Let r' be the graph obtained from T by gluing the vertices Vq and Vi together into a single vertex 
V, so that £y = £yg U Sy-^, and endowed with the 6-type condition 

(29) E^(^) = ("o + "i)/(^)- 
Then the eigenvalues of the two graphs satisfy the inequalities 

(30) A„(r)<A„(r')<A„+i(r). 

Proof. Similarly to the proof of Theorem 15.11 we consider the quadratic forms of the two graphs and 
observe that they are defined by exactly the same expression (see section . However, joining the 
vertices together imposes an additional restriction on the domain of the quadratic form of the graph 
r', namely 

D{h') = {f e D{h) : f{vo) = f{v^)}. 

Thus, the domain D{h') is a co-dimension one subspace of D{h) and the rest of the proof is identical 
to the proofs of the second and third inequalities in fl7r|) . □ 

Notice that if the domain of h' had co-dimension k, then one would have obtained by applying the 
same argument the inequality 

An(r) < A„(r) < A„+,(r). 

This observation immediately leads to the following generalization of Theorem 15.31 

Theorem 5.4. Let the graph T' be obtained from T by k identifications, for example by gluing vertices 
Vq, Vi, . . . f fc into one, or pairwise gluing of k pairs of vertices. Each identification results also in 
adding parameters aj in the vertex 6-type conditions, as in ^29\j. Then 

(31) A„(r) < A„(r) < A„+,(r). 

This statement can be also proved by the repeated application of Theorem 15.31 

Remark 5.5. Theorem 15.11 applied to the case a' = oo is also a result about joining vertices 
together, since the vertex Dirichlet condition has the effect of disconnecting the edges at the vertex. 
Note the difference with the result in Theorem 15.41 the eigenvalues of the joined graph are now 
shifted down, but not further than the next eigenvalue of F, which contrasts with the weaker "not 
further than fc-th next eigenvalue" result of the Theorem 15.41 

6. Some applications 

6.1. Dependence of the spectrum on the coupling constant a at one vertex. As in section [5] 
we consider the family of compact graphs Fq with a distinguished vertex v. Arbitrary self-adjoint 
conditions are fixed at all vertices other than f , while v is endowed with the 5-type condition with 
coefficient a: 

f is continuous at v and 
E £:{v) = af{v). 

We will be using the second condition in the form 

(32) (^ + l)E^(^)=i(^-l)/(^)' 

where \z\ = 1. The Dirichlet condition corresponds to z = —1, while the Neumann-Kirchhoff 
corresponds to z = 1. 

We denote by 'H{z) the operator with the above condition at v and the previously fixed set of 
conditions at all other vertices. Then main result of this section is representation of the spectrum of 
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this operator as the range at z of an irreducible multiple valued analytic function defined near the 
unit circle, plus a fixed discrete set, see Fig. [H 




Re"z 



Figure 1. Spectrum of a graph as a function of the parameter ^, which controls the 
matching condition of the form (152]) at one of the vertices of the graph. The figure 
shows the non-decreasing function A (2;) spiraling upwards and one of the constant 
branches from the set A. 



Theorem 6.1. There exist a bounded from below discrete set A G IR and a multiple valued function 
A{z) ( "dispersion relation") analytic in a neighborhood of the unit circle and real on , such 
that: 

(1) For any z G , one has 

(33) a(H(^))=A(z)UA, 

where oi^A) denotes the spectrum of the operator A. 

(2) The function A is irreducible, i.e. any of its branches determines by analytic continuation the 
whole function A. 

(3) Each analytic branch of A is monotonically increasing in the counter-clockwise direction of 
S\ 

Proof. Notice that, as we already know, for each z G 5*^, Hiz) is a bounded from below self-adjoint 
operator with a discrete spectrum. Let Xn{z) be its nth eigenvalue, counted with multiplicity in 
non- decreasing order. Then, according to the perturbation theory, it is continuous on the unit circle 
cut at z = —1, where one can observe that the ground state Ai has an one-side limit equal to —00. 
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Lemma 6.2. 

(1) The function A„(z) is non- decreasing in the counter-clockwise direction on the unit circle. 

(2) There exists a function F{X, z) analytic in a neighborhood of the unit circle such that 
A G a{l-i{z)) if and only if F{\, z) = 0. 

Indeed, the first claim of tlie Lemma is a replirasing of a part of Theorem 15.11 (and, essentially, a 
consequence of monotonicity of the quadratic form (fT6|) with respect to a). The second claim is a 
direct consequence of Theorem 13. 8[ 

Lemma 6.3. // either 

(1) A is a multiple eigenvalue of'H{z), for some z & S^, or 

(2) A belongs to the spectra of7i{zi) and 1-L{z2) for two different points zi,Z2 G S^, 
then A G aiTi^z)) for any z & S^. 

Indeed, if the eigenspace of 'H(z) corresponding to A is more than one-dimensional, then it contains 
an eigenfunction vanishing at the vertex v. Then, assuming z ^ —1, we conclude from ( l32l) that the 
sum of derivatives is zero as well. Therefore this eigenfunction works for any point z & S^. In the 
case z = —1 we select an eigenfunction with vanishing sum of derivatives at v and it automatically 
vanishes at v since it satisfies the Dirichlet condition. 

Let now A G criTi^zi)) f]o'(T-L{z2)) with the corresponding eigenfunctions /i and /2. Without loss 
of generality, zi ^ —1. First consider the case fi{v) ^ 0. Then, by Theorem 15. A as an eigenvalue of 
'K{z\) hes strictly between consecutive eigenvalues of 'K{z2) which contradicts our assumption. Thus 
/i(f) = and the sum of derivatives of f\ is zero by ( 132|) . We conclude that f\ is an eigenfunction 
for any z ^ thus completing the proof of the Lemma. 

We now define the set A as the set of all A G M such that A G cr('H(2;)) for all z ^ . As the 
previous lemma shows, one can describe A as the intersection (t('H(— 1)) ^(^{(l)) of the Dirichlet 
and Neumann spectra and thus is discrete and bounded from below. 

Consider now a neighborhood ?7 in C of the circle and the set S that is the closure of the set 



To put things simply, we remove from the union of the spectra all "horizontal" branches (z. A) with 
A G A. The second statement of Lemma 16.21 then implies that the set S is analytic. Moreover, for 
the neighborhood IJ being sufficiently small, it is a smooth complex analytic curve. Indeed, Lemma 
16.31 implies that for all A ^ A the eigenvalues are simple, and thus the eigenvalue branch is analytic. 
Let now A G A and (2:0, A) G S for some zq G . Then Lemma [6.31 says that the eigenvalue A near 
Zq either stays constant, or splits into the constant and possibly one or more increasing branches. 
The constant branches are excluded in the definition of S, equation Also, there can be no more 
than one increasing branch, otherwise one would find two values of z with the same values of A in 
E, which the lemma allows only to happen on horizontal branches. Moreover, as we have already 
concluded, the eigenvalue on such a branch must be simple. Then Rellich's theorem (e.g., [25] v.4. 
Theorem XII. 3]) says that the increasing branch is analytic. We thus conclude that S consists of one 
or more non-intersecting analytic curves. 

We will show now that there is only one component, if the neighborhood V is sufficiently small. 
First of all, the projection onto the A-axis of {\{S^ x M), where is any of the components of 
E is the whole real axis, otherwise the component whose projection does not cover the whole real 
axis would create an accumulation that would contradict the discreteness of the spectrum for each 
z G S"^. Then existence of two or more components would create equal eigenvalues for at least two 
distinct values of which would contradict to the Lemma 16.31 



(34) 
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Thus, S is an irreducible analytic curve that intersects each line {z} x M on the cylinder along the 
discrete spectrum o"('H(z)) and such that it is monotonically increasing counterclockwise. Hence, it 
forms a kind of a spiral winding around the cylinder infinitely many times when A — )■ oo and having 
a vertical asymptote when A — ?■ — oo (otherwise one would get a contradiction with the boundedness 
of each of the spectra from below). This proves the statement of the theorem. □ 

6.2. Nodal count on trees. Here we present a simple proof of a result that was first discovered 
in [Tl[23l[26] . With all ground-setting work done in Section [5] the proof is significantly shortened. 

Theorem 6.4. Let Xn be a simple eigenvalue of — -^ + V{x) on a tree graph T and its eigenfunction 
Z*^") be non-zero at all vertices ofT. Then f^"'^ has n — 1 zeros on T. 

Proof. To simplify notation we will drop the script n when talking about the eigen-pair (A„, Z*-""^). 
We will prove the result by induction on the number of internal vertices of the tree F. If there 
are no internal vertices, F is simply an interval and the statement reduces to the classical Sturm's 
Oscillation Theorem (see, e.g. [5]). 

Let f be a vertex of degree d^. Then v separates F into d^ sub-trees Fj, each with a strictly smaller 
number of internal vertices. For each subtree, the vertex f is a vertex of degree 1 with, so far, no 
vertex condition. We will impose 5-type condition with the parameter < oo chosen in such a way 
that the function / restricted to Fj is still an eigenfunction. This value is simply = fl{v) / fiiy). 
The eigenfunction / still corresponds to the eigenvalue A„ which is still simple (otherwise we can 
construct another eigenfunction for the entire tree F, contradicting the simplicity of A). It is an n^-th 
eigenvalue of Fj and, applying the inductive hypothesis, we conclude that the function / has — 1 
zeros on the subtree Fj. Thus we need to understand the relationship between the numbers nj and 
the number n. 

To do it, consider the subtree Fj oo which now has Dirichlet condition on the vertex v. By Theo- 
rem EUl there are exactly nj — 1 eigenvalues of Fj oo that are smaller than A. Now consider the tree 
Foo, which is the original tree but with the Dirichlet condition at the vertex v. On one hand it has 
exactly n — 1 eigenvalues that are smaller than A. On the other hand F^o is just a disjoint collection 
of subtrees Fj oo and its spectrum is the superposition of the spectra of Fj^oo- Therefore the number 
of eigenvalues that are smaller than A is 

n - 1 = ^(nj - 1). 
1=1 

Since we already proved that the number of zeros of / is equal to the sum on the right-hand side, 
we can conclude that yu(/) = n — 1. □ 
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Appendix A. Analytic Fredholm alternative 
Here we formulate the following version of the analytic Fredholm alternative: 
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Theorem A.l. Let A{z) be an analytic family of Fredholm operators of index zero acting in a 
Banach space E, where z runs over a holomorphically convex domain VL in C"" (or a more general 
Stein analytic manifold). Then there exists an analytic function f{z) in fl such that the set of all 
points z for which A{z) is not invertible coincides with the set of all zeros of the function f . 

The proof of this statement can be found in many places, e.g. it foUows from the Corollary to 
Theorem 4.11 in ^27\ . 
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